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0\ ■ Abstract 

We give the exact critical frontier of the Potts model on bowtie lattices. For the case of q = 1, the 
Q , critical frontier yields the thresholds of bond percolation on these lattices, which are exactly consistent 
with the results given by Ziff et al [J. Phys. A 39, 15083 (2006)]. For the g = 2 Potts model on the 
jrt bowtie-A lattice, the critical point is in agreement with that of the Ising model on this lattice, which 

has been exactly solved. Furthermore, we do extensive Monte Carlo simulations of Potts model on the 
bowtie-A lattice with noninteger q. Our numerical results, which are accurate up to 7 significant digits, are 



fi: 

i-rt ' consistent with the theoretical predictions. We also simulate the site percolation on the bowtie-A lattice, 

C^ - 

Q ■ and the threshold is Sc = 0.5479148(7). In the simulations of bond percolation and site percolation, we 



find that the shape-dependent properties of the percolation model on the bowtie-A lattice are somewhat 



^ ' different from those of an isotropic lattice, which may be caused by the anisotropy of the lattice. 
(N| ■ PACS numbers: 05.50.+q, 64.60. Cn, 64.60.De, 75.10.Hk 
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I. INTRODUCTION 

The reduced Hamiltonian of g-state Potts modelpj, 12| can be written as 

-n/kBT = KY,^'^.'^p (1) 

where ai, aj are the Potts spins on sites i and j. The Potts spin can take q values, namely 
(Tj = 0, 1, ■ ■ ■ , g — 1. The sum takes over all nearest-neighboring sites < i,j >. K is the coupling 
constant between spin ai and aj. This model can be mapped to the random-cluster model 3|, |4| 
with partition sum 



Q 



Y.v'^^q''^ (2) 



where ^ is a graph on lattice C that consists of Ni, random occupied bonds and Nc clusters. 
V = exp K — 1 is the bond weight, q is the cluster weight, which is not restricted to be integer 
in ([2]). An edge is set to be occupied by a 'bond' with probability p = v/{l + v) or vacant with 
probability 1 — p. A cluster is defined as a group of sites connected by the bonds. 

When q = 1, the random-cluster model ([2]) reduces to the bond percolation model [5j, in which 
the physicists have special interests. A central problem of the research in percolation model is the 
determination of the thresholds of the percolation transition on different lattices. After decades of 
investigation, the percolation thresholds on a lots of lattices have been determined. Most of the 
results are numerical ones, exact results are scarce. An important exact result for bond percolation 
threshold is the one for the triangular- type lattice (see Fig. [I^a)). As is summarized by ZifF et 
al. [g, lZJ , the percolation threshold of the lattice is determined by 



P{i,j,k)=P{iJ,k), (3) 

where P{i,j, k) is the probability that the three apexes i, j and k of the triangular cell (Fig. [11(b)) 
are connected in a same cluster via the bonds in the cell, while P{i,j, k) is the probability that the 
three apexes are isolated in three different clusters. There are other expressions to determine the 
percolation thresholds of the triangular-type lattices 6|], but they are in fact equivalent to ([3]). ([3]) 
was used to determine the exact percolation thresholds of a series of lattices, such as the martini 
lattices [6|,l8| and the stack-of-triangle lattices[9|]. All of these lattices belong to the triangular-type 
lattice, with different internal structure in the triangular cells. Ziff and Scullard 
to a series of lattices which they call bowtie lattices, as shown in Fig. [2 
percolation on the bowtie- A lattice was first determined by Wierman|ll 



transformation 17 1 . 



10] extend (|3]) 
The threshold of bond 
, using the star-triangle 




FIG. 1. (a) the triangular- type lattice, (b) the triangular cell, (c) the kagome-type lattice. The hatched 
triangular cells can have internal structures. 





(C) 




FIG. 2. Bowtie lattices: (a) bowtie-A lattice, (b) bowtie-B lattice, (c) bowtie-C lattice, (d) bowtie-D 
lattice. 



When g 7^ 1, the random-cluster model ([2]) can be considered as a general percolation model 
with cluster weight in the partition sum. The determination of the critical points of this model 
on different lattices is also a challenge in statistical physics. Similar to the problem of percolation 
threshold, exact critical points of the random-cluster model are also scarce. Up to now, exact 



critical points were found for only a small number o 



lattice, the honeycomb lattice, the martini lattices 12|, and the stack-of-triangle lattices 13|, [ij]. In 



two-dimensional lattices, such as the square 



these lattices, the honeycomb lattice and the martini lattices belong to the triangular-type lattice, 
with complex internal structure in the triangular cells. 

For the Potts model on the triangular-type lattice, the partition sum can be written as 



{a} A 



(4) 



where i,j, /c are the three apexes of the cell, and the sum and product take over all the cells. 



W^a(^, J, k) is the Boltzman weight of a hatched triangular cell, it can be written in the form of 

WAihJ, k) = A + B{5^^^^ + (5,^,,, + 5,,^J + C6^^^^,^. (5) 

For example, if the internal structure of the cells are simple triangles, then 

WA{i,j, k) = exp[K{6a,a, + S^.a^ + Sa^a,)] 

= (1 + v6^^^^){l + v6^^„J{l + v6^^„J 

= 1 + v{6^^^^ + 6^^^^ + 5.,^J + {v' + 3vy^^^^^^, (6) 

where v = exp K — 1 and 5o--o- o-^. = Sa-icr^cr a^- This gives A = 1,B = v, and C = f ^ + 3f ^. 



Kelland 



15| showed that the partition sum (jl]) is a self-dual one with self-dual point 

qA = a (7) 



Wu and Lin[l6|] established rigorously that this point is indeed the critical frontier in the ferro- 
magnetic regime 

35 + C > 0, 2B + C>0. (8) 

This is a very important result in determining the critical point of Potts model. For example, 
substituting the expressions A = 1 and C = v^ + 3v^ to ([7]), we get the critical frontier of the 
Potts model on the triangular lattice, which writes v^ + 3f^ = q. Substituting p = v/{l + v) to 
this equation with g = 1, we get the threshold of bond percolation on the triangular lattice, i.e.. 



Pc = 2sin(7r/18), which is a famous result [17 1. 



Basing on ([7]), Wu gives the critical points of the kagome-type lattices (Fig. Wic)) through a 
homogeneous conjecture |l3l| 

{q^A + 3qB + Cf - 3{qB + C)^ - {q - 2)C^ = 0. (9) 

This conjecture was checked by a lot of accurate numerical results given by the finite-size scaling 
analysis |14l|. which shows that Wu's conjecture was not exact but an excellent approximation. 

In this paper, we will show that (jTj) can be extend to the bowtie lattices for determining the exact 
critical points of the random-cluster model on these lattices. This is very similar to Ziff 's extension 
of ([3]) to this type of lattices for determining the bond percolation thresholds. Furthermore, we 
do Monte Carlo simulations of the random-cluster model on the bowtie-A lattice, the numerical 
results of the critical points are consistent with the theoretical predictions. In the procedure of 
simulations, we find that the shape-dependent properties [l8|, [l9| of the percolation model on the 
bowtie-A lattice are somewhat different from those of an isotropic lattice, by comparing some 
numerical results of bond percolation and site percolation. 



II. CRITICAL FRONTIER ON BOWTIE LATTICES 



In fact, ([3]) is exactly equivalent to the q = 1 case of ([7]). A is the sum of the weights of 
sub-graphs with the three apexes i,j, k are isolated in three different clusters, corresponding to the 
probability P{i,j, k). C is the sum of the weights of sub-graphs with the three apexes connected 
in a same cluster, corresponding to the probability P{i,j,k). For example, in a simple triangle, 
P{i,j,k) = {1— p)^ and P{i,j,k) = p^ + 3p'^{l—p). Substituting p = t>/(l+t>) to these expressions, 
one can easily find that P{i, j, k) and P{i, j, k) are different to A and C respectively with a constant 
[1/(1 + vyf . This property holds for any triangular cell (but the constant may be different for 
different internal structure). In Ref. [ij], a computer program was used to calculate the expressions 
of A, 5, C of n X n stack-of-triangle cell up to n = 4. In Ref. 9[, Haji-Akbari and Ziff give the 
expressions of P{i,j,k) and P{i,j,k). Comparing these expressions, we can see more clearly the 
equivalence of (E]) and the g = 1 case of d ?]). 

([3]) follows the duality arguments |20l. l2l| . it can be used to calculate the threshold of bond 
percolation on a lattice composed of triangular cells, and the lattice must be self-dual under the 
transformation shown in Fig. |3](a). The triangular- type lattice obviously meets this condition. 
Another type of lattice that meets this condition is the lattice shown in Fig. |3]^b), and the four 
bowtie lattices can be generated by inserting different internal structures in the hatched triangular 
cells of this lattice. This lattice is first given by Ziff and Scullard 10|], we call this lattice the 
'bowtie-generating lattice', because it 'generates' the four bowtie lattices. If simple triangles are 
inserted in the cells of the bowtie-generating lattice, it generates a lattice with double bonds as 
shown in Fig. HI The bond percolation with uniform probability p on the bowtie-A lattice can be 
got by setting the double bonds with probability 1 — a/1 — p, which is the Wierman's trick 111 ]. 





(a) 



(b) 



FIG. 3. (a) duality transformation of a triangular cell, (b) the bowtie-generating lattice, which is self-dual 
under the duality transformations of the cells. Four bowtie lattices as shown in Fig. [2] can be generated 
by this lattice by inserting different internal structures in the hatched triangular cells. 




FIG. 4. A lattice constructed by inserting simple triangles in the cells of bowtie-generating lattice and its 
bond weights of the random-cluster model (top). 

Similar to ([3]), ([7]) also follows the duality arguments [l5|, |l6|, this inspires us to extend ([7]) to 
calculate the critical points of random-cluster model on other lattice which meets the condition of 
self-duality. As an example, we calculate the critical points of the bowtie lattices in this paper. 



A. bowtie-A lattice 

For generality, we set the coupling constants of the three bonds of the triangle in the double- 
bond lattice with K, M, N. The bond weights are v = exp i^ — 1, u = exp M — 1 and w = exp A^ — 1, 
which are shown in Fig. |H The expressions of A and C are 

A = 1, (10) 

C = uvw + uv + vw + wu. (11) 

Substituting A and C to (jTj), we get the critical frontier of the Potts model on the double-bond 
lattice 

uvw + uv + vw + wu = q. (12) 

Setting 

M = K,N = K/2, (13) 

namely 



u = v,w = y/T+v — 1, (14) 
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we get a Potts model with uniform coupling constant K or random-cluster model with uniform 
bond weight v on the bowtie-A lattice. Substituting flT^ to ( IT^ . we get the critical frontier of the 
Potts model on the bowtie-A lattice, which writes 



(v"^ + 2v){\fTTv - 1) + w^ = g, bowtie-A lattice. (15) 

Solving (dS]) with g = 1, we get v^ = 0.679 312 786 or K^ = 0.518 384 653. This yields the threshold 
of bond percolation on the bowtie-A lattice, i.e., Pc = ^c/(l + 'Vc) = 1 —exp{—Kc) = 0.404 518 319, 



which is exactly equal to the result that given by Wierman 11| or Ziff et al. 10 |. 

The case of g = 2, which is the Ising model, has been exactly solved by Strecka and Canova[22| 
through mapping the model into a free-fermion eight- vertex model, and the critical points is given 
by 

[cosh(4is:/) - l]e^^' = 2, (16) 

where Kj = K/2 is the coupling constant of nearest-neighboring Ising spins. Solving f lT6|) . we get 
Kic = 0.333 135 959. Furthermore, if we solve ([15]) with g = 2, we get K^ = 0.666 271 918. The 
two results are exactly consistent with each other. 

The cases of g = 1 and 2 show the correctness of flTSj) for the bowtie-A lattice. For the other 
case, especially the case that g is not an integer, the correctness of ( 1T5|) is not verified. Thus we 
do extensive Monte Carlo simulation and finite-size scaling analysis for the Potts model on the 
bowtie-A lattice, the numerical results will be given in Sec. IIIII 

Setting the double-bond weights w = (namely A^ = 0), ( 1T2|) gives the critical frontier of Potts 
model on the square lattice 

uv = g, (17) 



which is a famous result 23 1 



B. bowtie-B, and bowtie-D lattices 

The critical frontier of the Potts model on the bowtie-B, the bowtie-C and the bowtie-D lattices 
can be obtained by the same way that we used for the bowtie-A lattice. In this paper, we only 
pay attention to the random-cluster model on these lattices with uniform bond weight v. 

In order to get the critical frontier of the Potts model on the bowtie-B and the bowtie-D 
lattices, we insert the network as shown in Fig. |5]^a) to the cells of the bowtie-generating lattice, 



y 



Vi 



\ 


t> 


V, ^ 


^ 


^ 


H^ 


L 


'v 


\ 




V3 


(a) 


(b) 



FIG. 5. (a) internal structure for constructing the bowtie-B and bowtie-D lattices, (b) internal structure 
for constructing the lattice as shown in Fig. [6^a). v = ex.p K — l,Vi = exp Ki — l,Wi = expLj — 1 are 
bond weights, M is the 3-site coupling constant between the three spins of the triangle in center. 

this generates a lattice that also has double bonds. The expression of A and C are 

A = q'^ + 5qv + 8v'^ + v^, (18) 



4 1 „.5\ 



C = 3v^ + v^ + vi{2qv' + 8v^ + 5v^ + v 



(19) 



where v = exp K — 1 and vi = exp Ki — 1 are the bond weights (see Fig. [5t^a)). Substituting them 
to (jTj), we get the critical frontier 



3v^ + v^ + Vi{2qv^ + SV^ + 5v^ + v^) - q{q^ + 5qv + 8v^ + w^) = 0. 



(20) 



Setting Ki = K/2, namely vi = a/1 + v — 1, we get the critical frontier of the Potts model on the 
bowtie-B lattice with uniform coupling constant K 



3v^ + v^ + {VTT^ - l)(2gf2 + 8v^ + 5v^ + v^) - q{q^ + hqv + Sv"^ + v^) = 0, 

bowtie-B lattice. (21) 

Solving this formula with g = 1, we get the critical point Vc = 1.142 305 296 or Kc = 0.761 882 490. 
This yields the threshold of bond percolation p^ = 0.533 213 122, which is exactly consistent with 
the value given in Ref. [lO(]. 

Setting K = 0, namely f 1 = 0, we get the critical frontier of the Potts model on the bowtie-D 
lattice with uniform K 



3v^ + v^ - q{q^ + ^qv + 8t;^ + v'^) = 0, bowtie-D lattice. 



(22) 



Solving this formula with g = 1, we get the critical point Vc = 1.669 919 079 or Kc = 0.982 048 164. 
This yields the bond percolation threshold pc = 0.625 456 813, which is exactly consistent with 
the value given in Ref. llOl . (I22l) is also the critical frontier of the Potts model on the martini-A 

n 

lattice, see Ref. [12[ for details. 
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C. bowtie-C lattice 





FIG. 6. After the decimation of the dotted spins, the lattice (a) becomes the lattice (b), namely the 
bowtie-C lattice. 



By inserting the network as shown in Fig. [5]^b) to the triangular cells of the bowtie-generating 
lattice, we get a lattice as shown in Fig . EJ^a). The expressions of A,C for the network as shown 
in Fig. [5]^b) have been given by Wu 12|: 



A = ■Uif2'U3 + ViV2{q + Wi+ W2) + V2V3{q + W2 + W3) + V3Vi{q + ^3 + Wi) 

+ {q + t;i + t;2 + fa) x [g^ + q(wi + W2 + W3) + h] (23) 

C = viV2V3h (24) 



where 



e^' -l,u;,- = e^" - 1 



h 



M+L1+L2+L3 



0L2 



- e^^ - e^^ - e^^ + 2 



L3 



(25) 



Ki, Li are the 2-site coupling constant, M is the 3-site coupling constant. Setting the bond weights 
wi = W2 = W3 = V2 = V3 = V (i.e., Li = L2 = L3 = K2 = K3 = K) and the 3-site coupling constant 
M = 0, we get the critical frontier of the Potts model on the lattice as shown in Fig. [6t^a): 



q{q + 2v)[q^ + 'iqv + t;^(4 + v)] + [q^ + hq^v - t;^(3 + v) + qv'^iS + v)]vi = 0. 



(26) 



Decimating the dotted spins in the lattice as shown in Fig. [6]^a), we get the lattice as shown in 
Fig. int^b), which is the bowtie-C lattice. If the bond weight vi is set as f 1 = w + a/w ^ + vq, we get 
the critical frontier of the Potts model on the bowtie-C lattice with uniform bond weight v. 



q{q + 2v)[q^ + 3gt; + f^(4 + v)] + [q^ + ^q^v - v^{'i + v) + gf^(8 + v)]{v + ^/v^^^+~qv) = 0. 

bowtie-C lattice. (27) 



Solving this formula with g = 1, we get the critical point v^ = 2.057 439 254 or Kc = 1.117 577 720. 
This yields the bond percolation threshold Pc = 0.672 928 906, which is exactly consistent with 
the result that given in Ref. |lOl |. 

We also solve (TT5|) . (pTI) . (l22ll . and (!27|) with other value of q, some results are listed in Table. [B 
In the table, we only give the results with g < 4, because the phase transition of the Potts model 
becomes discontinuous when q > 4 [23]. 

TABLE I. Exact critical points of the Potts model on the bowtie lattices, A = bowtie-A lattice, B = 
bowtie-B lattice, C = bowtie-C lattice, D = bowtie-D lattice. 



q 


Kc{A) 


Kc{B) 


Kc{C) 


Kc{B) 


1.0 


0.518 384 653 


0.761 882 490 


1.117 577 720 


0.982 048 164 


1.5 


0.601 796 427 


0.880 501 865 


1.263 365 870 


1.116 885 906 


2.0 


0.666 271 918 


0.971 526 940 


1.372 908 989 


1.218 755 726 


2.5 


0.719 275 592 


1.045 947 678 


1.461 235 663 


1.301 189 204 


3.0 


0.764 490 171 


1.109 155 997 


1.535 511 667 


1.370 686 799 


3.5 


0.804 030 933 


1.164 234 324 


1.599 748 425 


1.430 907 221 


4.0 


0.839 235 249 


1.213 123 525 


1.656 430 741 


1.484 126 262 



D. bowtie-dual lattices 



The critical frontier of the Potts model on the dual of the bowtie lattice, namely the bowtie-dual 
lattice (see Fig. 5 of Ref. [lO|), can be obtained by the dual relation 



m 



V ■ V 



Q, 



(2^ 



or 



-l)(e^ 



Q: 



(29) 



where K is the coupling constant on the bowtie lattice, and K* is the coupling constant on the 
dual lattice, v = exp K — 1 and v* = exp K* — 1 are the bond weights on the bowtie lattice and the 
bowtie-dual lattice, respectively. Thus the critical frontier of the Potts model on the bowtie-dual 
lattice is given by 



.K* 



l)(e-^-l) = g. 



(30) 
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III. MONTE CARLO SIMULATION AND NUMERICAL RESULTS 



A. Simulation of the Potts model on the bowtie-A lattice 



Using the cluster algorithm 2J| , we do Monte Carlo simulation of the Potts model on the bowtie- 
A lattice with a > 1. This cluster algorithm is a combination of Swendsen-Wang algorithm 25| and 



'coloring trick' 



26 



27(1, which drastically reduces the critical slowing down problem. This type of 



algorithm has been developed and used to simulate 0(n) loop model[28|, |29|, Eulerian bond-cubic 
model SOj and so forth. 

In the random-cluster model (|5]), the behaviors of the clusters are similar to those of percolation 
model. There are small clusters on the configurations when the bond weight v is small, but the 
clusters will growth as the value of v increases. The largest cluster diverges at the critical point Vc 
(or Kc) in the thermodynamic limit, which is called a 'percolating cluster'. However, for a finite 
system, the largest cluster may be different to the 'percolating cluster'. There are different rules 
to define the 'percolating cluster' for a finite system with periodic boundary condition. Here we 
use the 'wrapping cluster' SlJ, which is defined as a cluster that connects itself along at least one 
of the periodic directions. 

Basing on the definition of the 'wrapping cluster', the 'wrapping probabilities' on a rectangular 
lattice are defined as 



Rx = 


= (7^.), 


Ry - 


= (^.), 


? - 


= {n.y) 



(31) 

where (■ ■ ■ ) stands for ensemble average. T^a; is 1 (zero) if there is a (no) cluster that wraps along 
the X direction, whether or not the cluster wraps along y direction. The definition of TZy is similar 
to TZx- T^xy is defined as TZ^y = T^x ■ T^y, whose value is 1 if and only if there exists a cluster that 
wraps along x and y directions at the same time. 

The bowtie-A lattice that we simulated in the paper can be viewed as a 'rectangular lattice' 
with internal structure in the rectangular cells, see Fig. [3 This rectangular lattice is obviously 
not symmetry between x and y directions; therefore the values of Rx and Ry will be not the same. 



which will be shown concrete 
finite-size scaling behavior 132 



y later. The wrapping probabilities Rx, Ry and Rxy have the similar 



33| 



R = Ro + ai{K- K,)Ly' + a2{K - K,)L^y' + ■■■ + hiL^^ + 1)21^' + ■ ■ 



(32) 
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rectangular cell 




FIG. 7. Definition of wrapping probabilities on the bowtie-A lattice. 

where R=Rx, Ry or Rxy. Rq is the value of R at the critical point Kc. In present paper, we 
call it 'critical wrapping probability', yt is the thermal exponent of the Potts model, and yi the 



correction-to-scaling exponent, which has negative va 



ue. Oj and bi are unknown parameters. The 



value of yt can be derived by Coulomb gas method 3^ or conformal invariance 35] 

3 

with y/q = —2cos{TTg), 1/2 < g < I, 



(33) 
(34) 



where g is the coupling constant of the Coulomb gas. 

In order to illustrate our numerical procedure, we take the g = 1.5 Potts model as an example. 
The cluster algorithm easily allows us to do meaningful simulations for system with linear size 
up to L = 256. After equilibrating the system, 10^ samples were taken for each value of K for 
8 < L < 64, and 1.4 x 10^ samples for 128 < L < 256. Figure. |8]is an illustration of Ry versus K 
with different L for g = 1.5 Potts model. Using the Levenberg-Marquardt least squares algorithm. 




0.67 



0.6012 0.6015 0.6018 0.6021 0.6024 

K 

FIG. 8. Ry versus K for various system size of g = 1.5 Potts model. All the error bars are much smaller 
than the data points, the lines between points are added for illustration purpose. 



we fit the data according to the finite-size-scaling formula fl32|) . The fitting yields Roy = 0.7070(2), 
i?ox = 0.4363(2), /?ox2/ = 0.3780(1), K^ = 0.6017964(6) and yt = 0.887(1). In these results, we can 



12 



see that the critical point Kc is consistent with our theoretical prediction Kc = 0.60179643, and 
the thermal exponent is consistent with the Coulomb gas prediction yt = 0.8867. 

We also simulate the cases of g = 1, 2.5, 3, 3.5 and 4. All the numerical results and the theoretical 
predictions are summarized in Table. [Ill In the fitting of the data for g = 4 Potts model, logarithmic 
corrections 36l-l38l| should be included. Instead of (1321) . we fit the data according to 



R = Ro + ai{K - K,)Ly^{\ogLf^ + a2{K - K,f L^^^ {log Lf^^ + ■■■ + &i ^"f^"f^ ' ^' 



logL ' logL'^^^^ 

with yi < 0. However, we can see (from the table) that the results for g = 4 Potts model are 
relatively inaccurate. This is because of the limited system size. For an accurate fitting with 
logarithmic corrections, data with larger system sizes are necessary. 

TABLE II. Critical properties of the Potts model on the bowtie-A lattice, T=Theoretical results, 
N=Numerical results. The value of Kc (N) can be obtained by fitting the data of Ry, Rx, or Rxy, 
here we list the best estimation. 



Q 


^c(T) 


i^c(N) 


ytiT) 


ytm 


Roy 


Rqx 


Roxy 


1 


0.518 384 653 


0.5183847(7) 


0.7500 


0.750(1) 


0.6490(1) 


0.3838(1) 


0.3161(1) 


1.5 


0.601 796 427 


0.6017964(6) 


0.8867 


0.887(1) 


0.7070(2) 


0.4363(2) 


0.3780(1) 


2.5 


0.719 275 592 


0.7192756(11) 


1.1018 


1.101(2) 


0.7750(1) 


0.5151(2) 


0.4729(2) 


3 


0.764 490 171 


0.7644906(9) 


1.2000 


1.197(6) 


0.7975(1) 


0.550(1) 


0.515(1) 


3.5 


0.804 030 933 


0.8040314(25) 


1.3050 


1.296(12) 


0.8156(2) 


0.585(1) 


0.557(1) 


4 


0.839 235 249 


0.839235(2) 


1.5000 


1.50(7) 


0.829(5) 


0.62(1) 


0.60(1) 



B. Simulation of site percolation 



We also simulate the site percolation on the bowtie-A lattice, the finite-size scaling analysis 
gives Sc = 0.5479148(7), yt = 0.750(1), Roy = 0.6090(1), Rq^ = 0.4286(1) and Ro^y = 0.3349(1). 
The value of Sc coincides with 0.5475(8), given by van der MarckJ3^], and our result is much more 
accurate. 



IV. SHAPE-DEPENDENT PROPERTIES OF THE BOWTIE-A LATTICE 



We also do Monte Carlo simulations of the bond and site percolation models on the square 
lattice. For the percolation thresholds, the numerical estimations give pc = 0.5000000(4) (bond), 
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and Sc = 0.5927460(5) (site). The accuracy of Sc is somewhat lower than the resuh given by 
Feng et al 40| . The numerical results of the critical wrapping probabilities are listed in Table IIIH 
these results give i?oe=-Rox + -Roj/ — -Roxj/=0. 6902(5) for the bond percolation and 0.6905(5) for the 



19ll4l|. Roe is 



site percolation, which are consistent with the literature value i?oe = 0.6904737 

the critical value of the probability that there is a cluster wrapping along one or both of the two 

coordinate directions. 




(a) 



(b) 



FIG. 9. (a) A rectangular lattice (on macroscopic level) with square symmetry on microscopic level, 
with aspect ratio r = 3/4; (b) A rectangular lattice (on macroscopic level) with triangular symmetry on 
microscopic level, with aspect ratio r = ^3/2. 



For a lattice that is isotropic on microscopic level (equivalent in the coordinate directions, or with 



the triangular or hexagonal symmetry, etc.), the va 



of the system shape and boundary conditions 
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ues of the wrapping probabilities are functions 



4l|, which are independent of the percolation 



type. In saying the shape, it is on the macroscopic level. For example, a macroscopic rectangular 
lattice may be a lattice with square or triangular symmetry on microscopic level, see Fig. |9] for 



examples or see Ref. 
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and 



19l for details. From now on, in saying 'isotropic' or 'anisotropic' in 
the text, it is on the microscopic level. For an isotropic rectangular lattice with periodic boundary 
conditions (not twisted), the value of the wrapping probability -Roe is determined by the aspect 
ratio of the lattice 19 1 

1 



R^Ar) 



-[Z,(8/3,r)-Z,(2/3,r)], 



(36) 



where r is the aspect ratio and Zc{h,r) is 



Zcih) 



rf{w) 



E E expj-^K^ + nV]} 



(37) 



n=— oo n'=— oo 



rj{w) 



,1/24 



k=l 

Roe = 0.6904737 for the square lattice. 



is the Dedekind eta function and w = e ^'^''". When r = 1, (1361) gives 
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In our Monte Carlo simulations, the bowtie-A lattice is viewed as a rectangular lattice (with 
periodic b.c), as shown in Fig. [TJ and the aspect ratio is r = 1/v^- However, the wrapping 
probability of the bowtie-A lattice is not determined by this aspect ratio. (1361) gives i?oe(l/v3) = 
0.7439918, this is inconsistent with our numerical values i?oe = 0.7167(1) (bond) and 0.7027(1) 
(site). Furthermore, we can see that the values of Rq (include Roe,Rox,Roy, and Roxy, see Table. 
mil ) for the bond percolation are different from those for the site percolation. These results show 
the difference between the bowtie-A lattice and an isotropic lattice in the aspect of shape-dependent 
properties. 

The difference may be caused by the anisotropy of the bowtie-A lattice, and ( l36l) is only valid 
for an isotropic lattice. After mapping the percolation model on the bowtie-A lattice onto the 
Gaussian model[4l|, it will be an anisotropic one. Thus a rescaling on x or y direction is required 
in order to obtain an isotropic Gaussian model. After the rescaling, the isotropic model would 
have an aspect ratio that is not equal to l/-\/3 (Unfortunately, we don't know the value of the 
effective aspect ratio). This is the reason why our numerical value of i?oe for the bowtie-A lattice 
doesn't coincide with _Roe(l/"\/3)- 

More important, there is no reason why the bond and site percolation models on the bowtie-A 
lattice should map onto the same (anisotropic) Gaussian model, and after rescaling to the isotropic 
Gaussian model, there would be different aspect ratios for the bond and site percolation models. 
Thus, there would be different values of the wrapping probabilities for the bond and site percolation 
models, which are verified by our numerical results. 

Another quantity that we simulate for investigating the shape-dependent properties is the av- 
erage density of clusters, which is defined as the average number of clusters per site. We find that 



on such an anisotropic lattice, the average density of clusters also behaves as|18| 

n = nc + b/N + --- , (38) 

when the system is at the critical point, ric is the value of n for the infinite system, and A^ the 
total number of sites of a finite system, b is the excess number of clusters over the bulk value of 
clusters ricN. For an isotropic lattice with periodic boundary conditions (not twisted), the value 
of b is also a function of system shape, which has the same value for the bond and site percolation 



models|l8l| 



b(r) = ^ + w'>/\2V3r - ^) + ^^(^3^ _ 1) + ^5/48 ^ 2u;53/48 _ ^23/16 ^ ^77/48 ^ . . . ^(gg) 

where w = e~^'^^. This equation gives b{l/\/3) = 0.943713, which is obviously inconsistent with 
our numerical results for the bowtie-A lattice: b^ = 0.9120(6) and b^ = 0.8957(4) (B denotes bond 
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percolation, S denotes site percolation). More important, the value of b^ is different from that of 
b^. The physical essence of these results for b is the same to that for the wrapping probabilities, it 
may also be caused by the anisotropy of the bowtie-A lattice. 

In conclusion (Table lllip . the shape- dependent properties of the percolation model on the 
bowtie-A lattice are somewhat different from those of an isotropic lattice, this may be caused 
by the anisotropy of the bowtie-A lattice. The bond and site percolation models on the bowtie- 
A lattice correspond with different anisotropies after mapping onto Gaussian models, and thus 
with different aspect ratios after rescaling to isotropic models, which lead to the difference of 
shape- dependent properties between the bond and site percolation models. 



TABLE III. Critical values of the wrapping probability Rq, the average density of clusters ric and the 
excess number of clusters b on the square and bowtie-A lattices, Site=Site percolation, Bond=Bond 
percolation. '^=Reference[l8], =Reference[19'],*^=Reference[4ll],'' =Reference[44] 



System 


Roy 


Rox 


Roxy 


Roe 


b 


ric 


Bond (square) 


0.5208(5) 


0.5210(5) 


0.3516(5) 


0.6902(5) 


0.883(1) 


0.098076(1) 




0.5210583'^'^ 


0.5210583'^'^ 


0.3516429^^'^ 


0.6904737"^'^ 


0.8838(5)^ 
0.883756 ^ 


0.0980763(8)^= 

3V3-5c 
2 


Site (square) 


0.5211(5)(1) 


0.5209(5) 


0.3514(4) 


0.6905(5) 


0.8834(11) 


0.027598(1) 




0.5210583'^'^ 


0.5210583'^'" 


0.3516429"'''= 


0.6904737"^'^ 


0.8832(3)^ 


0.0275981(3)'= 


Bond (bowtie-A) 


0.6490(1) 


0.3838(1) 


0.3161(1) 


0.7167(1) 


0.9120(6) 


0.119212(1) 


Site (bowtie-A) 


0.6090(1) 


0.4286(1) 


0.3349(1) 


0.7027(1) 


0.8957(4) 


0.023990(1) 



V. CONCLUSION 

In conclusion, we have given the exact critical frontier of the Potts model on the bowtie lattices. 
For q = 1, the critical frontier yields the thresholds of bond percolation on these lattice, which 
are consistent with the results given by Ziff et al. For the q = 2 Potts model on the bowtie-A 
lattice, the critical point coincides with the exactly solved critical point of the Ising model on this 
lattice. Furthermore, we do Monte Carlo simulations of the Potts model on the bowtie-A lattice, 
the numerical results, which includes the cases with noninteger q, are in agreement with the exact 
critical points in a high precision. 

Furthermore, by comparing the numerical results of the critical wrapping probability i?o smd 
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the excess number of clusters b for the bond and site percolation models, we find that the shape- 
dependent properties of the percolation model on the bowtie-A lattice are somewhat different from 
those of an isotropic lattice, which may be caused by the anisotropy of the bowtie-A lattice. 
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